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Abstract. We present a detailed study of the circumstellar envelope of the post-asymptotic giant branch "21 fim 
object" HD 56126. We build a detailed dust radiative transfer model of the circumstellar envelope in order to 
derive the dust composition and mass, and the mass-loss history of the star. To model the emission of the dust we 
use amorphous carbon, hydrogenated amorphous carbon, magnesium sulfide and titanium carbide. We present 
a detailed parametrisation of the optical properties of hydrogenated amorphous carbon as a function of H/C 
content. The mid-infrared imaging and spectroscopy is best reproduced by a single dust shell from 1.2 to 2.6" 
radius around the central star. This shell originates from a short period during which the mass-loss rate exceeded 
10"* M0/yr. We find that the strength of the "21" fim feature poses a problem for the TiC identification. The 
low abundance of Ti requires very high absorption cross-sections in the ultraviolet and visible wavelength range 
to explain the strength of the feature. Other nano-crystalline metal carbides should be considered as well. We 
find that hydrogenated amorphous carbon in radiative equilibrium with the local radiation field does not reach a 
high enough temperature to explain the strength of the 3.3—3.4 and 6—9 ^m hydrocarbon features relative to the 
11 — 17 fim hydrocarbon features. We propose that the carriers of these hydrocarbon features are not in radiative 
equilibrium but are transiently heated to high temperature. We find that 2 per cent of the dust mass is required 
to explain the strength of the "30" /im feature, which fits well within the measured atmospheric abundance of 
Mg and S. This further strengthens the MgS identification of the "30" fim feature. 

Key words. Stars: individual: HD 56126 - Stars: AGB and post-AGB - Stars: carbon - Circumstellar matter - 
Stars: mass-loss - Infrared: stars 



1. Introduction most recent evolutionary descendants of the AGB stars. 

■ ■ ' Mass loss near the tip of the AGB is not well understood. 

Post-asymptotic giant branch stars (post-AGBs) are ob- ^j^^^^ indications that the mass-loss rates near the 

jects in which the strong AGB mass-loss has ceased while ^j^^ ^GB exceed the values predicted by the exist- 

the remammg star is not (yet) hot enough to lomse the -^^ ^^^^ ^^.-^^^ ^-^^ ^^^^^^ CRL 2688). There is also 

surroundmg material. These objects are characterised by a increasing evidence that non-spherical mass loss is com- 

double-peaked spectral energy distribution (SED) . In the ^^^^ ^j^^ ^1^^ pj^^^^ g^^^^^ evolution 

ultraviolet (UV) and visible range the central star is visi- tij^es^ales are determined by mass loss, these poorly un- 

ble. In the mid-infrared (IR) range the light is dominated ^erstood phenomena limit our knowledge of the evolution 

by emission from the dusty circumstellar envelope (CSE) . ^^^^^ ^^.^^^^ subsequent stages. 

Post-AGB stars are the prime objects to study stellar evo- ^^^^ j^^^ determinations are directly coupled to the 

lution near the very end of the AGB because they are the composition of the newly condensed dust. The dust com- 

position of the CSE reflects the physical conditions in the 



outer atmospheres of the AGB star during the phase when 



Send offprint requests to: S.Hony, e-mail: shony@rssd.esa.int 
* based on observations taken at the European Southern 

Observatory, La SiUa, Chile and observation obtained with ^^^^ dust was formed and the processing that took place af- 

ISO, an ESA project with instruments funded by ESA Member terwards. We study the composition of the dust surround- 

states (especially the PI countries: France, Germany, the ing the post-AGB star HD 56126. This star is member of 

Netherlands and the United Kingdom) with the participation a class of rare objects that exhibit an emission featu re at 

of ISAS and NASA 21 ^m ((Kwok et al.iil98ai . The 12 fsee iKwok et al.iil999l 
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and references therein) post-AGB stars that show this fea- 
ture are commonly referred to as the "21 /xm objects" . The 
"21 fiia objects" form a group with rather homogeneous 
properti es: carbon-rich, low metallicity and s-process en- 
hanced llVan Winckel fc Revniersll20o'o|) with lar ge IR ex- 



cesse s ("Kw ok et al spectral types G or F (|Hrivna]3 

ll995iViVolk et alJ l)l999l) have shown that the shape of the 
"21" fim feature between different sources shows hardly 
any variation. Recently we have reported the dete ction 
of the sa me feature in two PNe ((Honv et al.l and 
IVolk et al. (2002) have reported on its detection in the 
PN IC 418. 

Besides the "21" fim feature the "21 fj,m objects" ex- 
hibit other infrared emission features that are unique to 
these environments. The carriers of many of these features 
have not been firmly identified and there is considerable 
debate in the literature on the composition of the dust 
around the "21 /xm objects". In the next paragraphs we 
shortly review these features and the proposed carriers. 

The carrier of the "21" fim feature has re- 
mained unidentified for more than a decade, al- 



OCfNH?), JSourisseau et al.lll992ft. SiS^ ( 


Goebell 


lBegemamietalJ'l9961. fuUeranes 1 Webster 


1995|). nano- 



bon (HAG) (IGrishko et alJl200lh and oxygen bearing side 
groups in coal l|Papoularl200CI|) . However, none of the spec- 
troscopic comparisons between the laboratory spectra of 
the suggested materials and the emission from the post- 
AGB stars is satisfactory. Recently we have suggested ti- 
tanium carbide nano- crystals (nano-TiC) as t he carrier of 
the "21" ^m feature l|von Helden et al.ll20n(]l) . This iden- 
tification is based on i) the excellent spectroscopic match 
of nano-TiC with the 21 /xm feature, ii) the fact that 
the nano-TiC resonance does not depend on the crys- 
tal size for the size-range for which measurements are 
available - in good agreement with its observed constant 
profile, in) the fact that TiC is found in presolar grains 
ijBernatowicz et al]ll996l) for which the isotopic ratios im- 
ply an origin in the C SE of carbon-rich evolved stars 
( Anders fc Zinner 



We note that the nano-TiC identification is not com- 
monly accepted either. The main counter argument is that 
the conditions (high density and pressure during dust for- 
mation) required to form TiC in the eje cta of carbon-rich 
evolved stars may not be met ijKwok et al...,2002'l . Another 
argument against nano-TiC is the strength of the "21" /im 
feature. Since Ti is a rare element, it may be difficult for 
TiC to produce a feature that stands out as much as it 
does. 

As explained above the "21 /xm objects" exhibit other 
emission features for which the carriers are still under de- 
bate. Besides the well-known unidentified infrared (UIR) 
bands at 3.3, 6.2, "7.7" and 11.2, generally attributed 
to polycyclic aromatic hydrocarbon (PAH) molecules 
ijAllamandola et al]ll989l) . these sources exhibit a 3.4 /im 
feature and broad emission features (so-called plateau 
features) between 6—9 and 11—17 /xm (^Buss et al.i. 199(1 



ljusttanont et alJll996j) . As a class, the emission bands ob- 
served from the "21 nm objects" are distinctly different 
from the UIR feat ures from most sources (see Sect. I3.l(l . 
iKwok et al] l|200lh have discussed these features in terms 
of alkane and alkene side groups on very l arge aromatic 
molecules and small carbonaceous particles. IGuillois et alJ 
l|l996^ have compared them to those measured in coal 
(a natural terrestrial carbonaceous compound). The spec- 
troscopic match with coal is very encouraging. However, 
a fundamental questions concerning the identification is 
whether carbonaceous grains can reach the temperatures 
required to emit strongly in the 6—9 /im range. 

Lastly, these sour ces exhibit a strong broad feature 
between 23— 45 Atm |[Omont et all Il995t ISzczerba et alJ 
i997; Hrivna k et all 1200 0), called the "30" /um feature. 
The carrier of the "30" /xm feature is identifie d with mag- 
nesium sulfide (MgS) fSee lHonv et alJ 120021 and refer- 
ences therein). The identification of the "30" /Ltm feature 
with MgS seems solid but open questions remain on the 
heating of the MgS grains and the amount of MgS required 
to explain the prominence of the feature. 

A meaningful approach towards the identification of 
the various dust components in the CSE of these post- 
AGBs is that of radiative transfer modelling of the SED 
of these objects. This may also help to improve our un- 
derstanding of the mass loss at the very end of the AGB 
phase. Such modelling allows a comparison of the emis- 
sion characteristics of candidate materials and CSE emis- 
sion and yields insight into the relative contributions of 
the different materials present in the CSE. 

We have elected to model HD 56126 (IRAS 
07134-1-1005, SAO 96709, BB+IO 1470) in detail. 
HD 56126 has the strongest "21" fim feature relative to 
the other dust features of all known "21 /im objects". 
Therefore, HD 56126 provides the strongest constraint on 
the abundances needed to explain the "21" /xm feature. 
This is even m ore so because the photosphere of the star 
is metal poor l|Van Winckel fc R.evnier I hnod . Both the 



central star and the dust shell of this source have been 
well studied in the literature. Especially important in this 
respect are images at various wavelengths i n the niid-IR 
in which the source is well resolved (M cixner et alJ|l997t 
iDaval et ablll 998l:I.Tura et"ini200ntrKwok et al. 1120021 and 
this work). These images are essential to constrain the 
model parameters. 

This paper is organised as follows: In Sect. [21 we give 
the stellar parameters of HD 56126. In Sect.|31we present 
the infrared measurements that we use. In Sect. ^ we con- 
sider the energy balance in the nebula and constrain the 
stellar effective temperature. In Sect. we list the param- 
eters that we use to model the absorption and emission in 
the CSE of HD 56126. The optical properties of the dust 
constituents are discussed in Sect. El In Sect.|Z|we present 
the model results. In Sect. (HI we discuss our findings and 
finally in Sect. El we give a summary and conclusions. 
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Table 1. UV to IR photometric data of HD 56126. Units Table 2. Details of the ISO observations used in this 



are magnitudes unless given explicitly. 



study. 



Walraven 

V V-B B-U U-W B-L 

-0.548 0.365 0.75 0.556 0.312 fPel 1976) 

-0.555 0.400 0.766 0.588 0.322 (vG 1986)" 
Johnson ^M^Gen^erei^^l^^S^ 
B V 
9.13 8.23 



Near-IR (Hrivnak et al. 1989) 
J H K L M 

6.92 6.66 6.65 6.43 6.10 



Stromgren 
b-y (v-b)-{b-y) {u-v)-{v-b) 
0.639 0.182 1.477 

[ergs s~^cm~^A~"'] 



fOlsen 19931 



u y 
1.03e-12 1.71e-12 



fUeta et al. 2000) 



Geneva (van Winckel, priv.comm.) 
U B Bl Ba Vi G 

11.010 8.337 9.452 9.616 8.970 9.213 



"van Genderen et al. (1986 



2. Stellar parameters 

HD 56126 is a bright star with a visual magnitude of 
8.3. The star has b een classified a s a F5 supe rgiant by 
iNassau et all l)l965j) . More recentlv iHrivnak et al. (198^ 
published a medium resolution optical spectrum and 
determined the spectral type to be F0-5I, consistent 
with the previous spectral classification. There are sev- 
eral photosp heric abundance determina t ions published for 
HP 56126 llParth asarathv et all Il993 iKlochkoval 11995: 
IVan Winckel fc R ev niers 2000 1). These studies all demon- 
strate the low photospheric metallicity of the star, with a 
value of [Fe/H] = —1.0 and enrichment of C, N, O and 
the s-process elements. There is less agreement on the 
exact amount of enr i chmen t of the CNO-elements. While 
IParthas arathv et alJ l)l992j) and lVan W inckel fc Revniers' 
(,2000 ) find a C/0 value close to unity, Klochkova (,1995) 
determine a lower abundance of O and derive a C/0 value 
of ~1.5. 

The conversion from spectral type to effective temper- 
ature (Toff) is not clear-cut. Due to the peculiar abun- 
dances one should be critical of applying standard spec- 
tral type calibrations based on solar type abundances. For 
example, the tables o fEZ Jager & Nieuwenhiiiize 
quote 6370 for a F5I star with solar abundances. Values of 
Toff for HD 56126 within a large range (5900-7250 K) are 
derived by various methods. We give most weight to the 
temperatures determined in detailed abundanc e analyses 
using high resolution spectroscopy (|Parthasar9 .1 ^hv et al.l 
ll992HVan Winckel fc R,evniersll200(Tl) . These studies yield 
a high value of the effective temperature (7000—7250 K). 
This temperature corresponds to a F1-2I spectral classi- 
ficatio n, still within the range derived by iHrivnak et alJ 
l|l989l) . Such a high value for Tcff also fits better with the 
observed spectral energy distribution (SED) (see below). 





HD 56126 ISO observations 


Instrument 


Obs." 


Q 5 TDT*" 




Mode 


(J2000) (J2000) 


sws 


06 


071610.20 +09 59 48.01 71802201 


sws 


06 


071610.30 -F09 59 48.01 72201702 


sws 


01(3) 


071610.20 +09 59 48.01 72201901 


LWS 


01 


071610.20 +09 59 47.80 72201802 



"Observing mode used see jde Graauw ct al. (119961) and 
IClegg et al.l lll996l) . Numbers in brackets correspond to the 
scanning speed. 

''TDT number which uniquely identifies each ISO observation. 



ljura et all l)200Cl|) have made the case that the low 
iron abundance of the star, the high galactic latitude 
(&==+10°) and its high radial velocity make it probable 
that HD 56126 is a population II star with a zero age 
main-sequence mass of ~1.1Mq. 

The distance {d) to HD 56126 is poorly k nown. Values 
betwe en 2 and 3 kpc are used in the literature. I Yuasa et alJ 
l|l999|) derive a statistical distance of 2.12 kpc from a 
principle component analysis based on IR AS photome- 
try an d the molecular expansion velocity. iKnapp et alJ 
l|200(]|) derive 2.4 kpc from the radial velocity assuming 
that the rad.ial ve locity is due to the galactic rotation while 
I.T^ira et al\ l(2000l) derive 2.3 kpc from the theoretical lumi- 
nosity of 6600 L© and the observed flux. In the following 
we elected to use (i=2.4 kpc. Note that this distance yields 
a luminosity of 6000 L(n which is higher than expec ted 
(-3000-4000 Lq) (e.g. lBoothrovd fc Sackma^l 19881) for 
a star with a ~1.1 M© ZAMS mass on the basis of the 
initial mass-final mass relation and core-mass-luminosity 
relation. In view of the uncertain distance, we explicitly 
give the distance dependence of the values that we derive. 



3. Observations and data reduction 

3.1. Infrared spectroscopy 

For the spectral energy distribution (SED) of HD 56126 
we use IR data obtaine d with the Shor t Wavelength 
Spectrometer (SWS) l|de Graauw et alJ 1 19961) and 
Long Wavelength Spectrometer (LWS) ijClegg et alJ 
Il996|) on-board th e Inf rared Space Observatory 
(ISO) l)Kessler et al.l Il996|) . These data consist of 
one SWS/AOTOl spectrum from 2.3-45 /im, two 
SWS/AOT06 spectra covering the ranges 16.5—24 and 
23—43 fim and a LWS spectrum from 45 to 200 /im. 
Details of the ISO observations are given in Table |21 

The SWS data reduction steps involved the removal of 
cosmic ray hits. Next, we combine all available data from 
the SWS observations in order to maximise the signal-to- 
noise ratio. We correct each separate observation, subband 
and detector combination in such a way that its average 
equals the average of all available data at the correspond- 
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Wavelength [|xm] 

Fig. 1. The spectral energy distribution of HD 56126. We show the lUE data (crosses), optical and near-IR photometry 
(triangles) (see Tabled for details), IRAS photon ietry (diamo i ids) a nd t he ISO/SWS and ISQ/L WS spectra (soHd 
line). In grey we show near-IR spectroscopy from IK wok et al.1 l)l99Cl() and IPudm aiier et al. For reference we 

also show a Kurucz model (dotted line) of the central star with T sfF=7250 K, logq=1.0 a nd [M/H] = -1.0 ijKurucd 
Il993l) reddened with E(B-V)=0.33 using a standard extinction law ijSavaee fc Matliislll979() . The standard extinction 
law is clearly not applicable. In grey we indicate the contribution of the "21" /zm feature and the dashed area indicates 
the "30" /Ltm feature due to magnesium sulfide. Note that on this scale surface corresponds to energy. Therefore one 
may readily see that an appreciable amount of energy is emitted in the "21" and "30" /zm features. The inset shows 
the Kurucz model before reddening (solid line) and the reddened model atmosphere (dotted line). 



ing wavelength. Below 12 /im we apply an offset and above 
12 ij,m a scaling to the data. We apply an additional scal- 
ing factor of 1.1 to the data of sub-band 3C (16—19.5 /im) 
to form a continuous spectrum. Finally, the data are re- 
binned to a fixed resolution grid of A/AA=200. 

Data reduction of the LWS observation consisted of 
extensive bad data removal, rebinning on a fixed resolu- 
tion grid of A/AA=50 and spHcing of the data to form 
a continuous spectrum. We apply scaling factors to the 
data from 45—70 /im and offsets at longer wavelengths. 
The data of detector 7 (130—140 /im) are very bad and 
we have replaced those by a constant level of 10 Jy, the 
average of the much better data on either side of detector 
7. The combined SWS /LWS spectrum together with data 
from the literature is presented in Fig. 

The SED of HD 56126 is double peaked. At wave- 
lengths shorter than 4 /im we see the stellar photosphere 
directly. At longer wavelengths the dust emission domi- 
nates. We see continuum emission which peaks at ~25 /im 
with strong broad emission features perched on top. In 



Fig. n we indicated the broad continuum due to amor- 
phous carbon (a-Carbon) and the "21" /xm, "30" /xm and 
the plateau feature, (tentatively) identified with titanium 
carbide, magnesium sulfide, and HAG or PAHs, respec- 
tively. 

In Fig. El we show a detailed view of the plateau fea- 
tures. We indicate the emission that can be attributed to 
aromatic and aliphatic bonds separately. Notice the dif- 
ferences with the emission spectrum of the planetary neb- 
ula NGC 7027 in feature shape and relative strength of 
the bands. The aliphatic component is much stronger in 
HD 56126 than in NGC 7027. 

3.2. N-band imaging 

We also obtai ned an 11.9 /im image with the TIMMI2 
l|Reimann et al. 2000: Kaufl et al. 2000 ) mid-IR camera 
attached to the ESO 3.6m telescope at La Silla, Chile on 
25 December 2001. The camera is equipped with a 320x240 
pixel array; we applied a pixel scale of 0.3"/pixel. We used 
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Fig. 2. The 3 to 17 /im spectrum of HD 56126 after re- 
moving the contribution of the star and the dust contin- 
uum. Several emission features due to vibrationally ex- 
cited C-H and C-C bonds in hydrogen-carbon compounds 
are identified. We have indicated the features due to aro- 
matic (grey lines) and aliphatic (black lines) structure. 
The inset shows a blow-up of the 3 /zm region, where the 
dominant aromatic 3.3 /im and the weaker 3.42 /^m signa- 
ture are observed. For reference we show the spectrum of 
the PN NGC 7027 as a prototypical PAH spectrum (grey). 
The emission features of HD 56126 are markedly different 
from the PAH spectrum of NGC 7027. 



a chop throw of 15" north-south and a nod of 15" east- 
west. This allowed for both the chopped as well as the nod- 
ded positions to fall onto the detector. We combined the 
positive and negative images resulting from the chopping 
and nodding positions. The images were reduced using 
a shift-and-add technique, where the shift between indi- 
vidual frames was determined from a least-squares com- 
parison between the individual frames. The resulting im- 
age was then deconvolved using an empirically determined 
point spread function obtained by observing HD 32887 
with the same set-up as HD 56126. We show the final 
image in Fig. |31 

The TIMMI2 11.9 /im image shows a clearly resolved 
envelope of ^^5" diameter. The emission peaks at 1.1" 
to the east from the centre of the nebulae. There are 
two other local intensity maxima: in the west and in the 
north, although the latter is not clearly resolved. Our 
11.9 /i m image is consistent with other mid- IR imaging 
studies (Meixner et al.'l997':'Dava l et alJll99llJura et alJ 
r2000t .Kwok et al. 2002) that show the same morphology 
throughout the 8—21 ^m wavelength range. The emis- 
sion maximu m in the north is resolved in the images of 
iKwok et al.l ()2002,'l . We derive the azimuthal averaged in- 




N 



11.9 |im 

Fig. 3. TIMMI2 image of HD 56126 in the 11.9 ^m narrow 
band filter. The scale on the top left indicates 1". There 
are three emission maxima on a ring-like structure (indi- 
cated by a dotted line) at ~1.1" from the centre (indicated 
by the '-I- '-sign). The grey-scale shows the square of the 
measured intensity to enhance the contrast of the emission 
maxima. The contours indicate the 5—95 per cent inten- 
sity levels in steps of 10 per cent. The outer contour (5 
per cent) lies at ~2.4" from the centre. 



tensity profile from Fig.|21to compare with the model out- 
put. 

4. Basic considerations 

Before we describe the details of our model we consider 
a few basic properties and the simple constraints that we 
can derive from them. 

4.1. Energy balance 

In Fig. ^ equal areas under the curve correspond to equal 
energy and it is clear that most of the energy is emitted 
in the IR. In Table 01 we specify the fractions of the en- 
ergy emitted by each of the components in Fig. ^ Only 
44 per cent of the light comes from the star directly while 
56 per cent has been absorbed in the CSE and re-emitted 
in the IR. This implies that the extinction averaged over 
all angles around the star must be appreciable. Using a 
standard extinction law and a 7250 K Kurucz model at- 
mosphere (Kurucz 1993) we estimate that the average Ay 
equals 1.1, which corresponds to rv=1.2. However, such 
a value for the visual extinction corresponds to hardly 
any extinction in the near and mid-IR. Thus the near- 
IR (2—4 fim) observations give the true continuum level 
of the star since no obvious dust contribution is detected 
in this wavelength range. Therefore, for a given effective 
temperature of the stellar photosphere there is no freedom 
in choosing the absolute scale. With this we can directly 
compare the integrated flux in the observed SED to the 
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Table 3. The integrated contributions of the differ- F8 F5 F2 FO 

ent components of the spectral energy distribution (see . ^ ^ ^ ^ 




integrated flux from the model atmosphere. In Fig. ^ we 
show the integrated flux in the complete observed SED 
(•^obs) divided by the integrated flux in a Kurucz model 
atmosphere (^atm) as a function of Tcff . As parameters of 
the Kurucz models we choose log 5= 1.0 and a metallicity 
of [Fe/H]=-1.0 as determined for this star. For Toff in the 
range 7000—7500 K ^obs/-^atm is approximately one. This 
ratio should be close to unity if the following assumptions 
hold, i) The extinction in the near-IR is negligible, ii) The 
extinction along our line of sight is similar to the average 
value of the extinction. 

The first assumption holds unless there is a large 
amount of grey extinction. Extinction caused by particles 
smaller than the wavelength at which they absorb falls 
off steeply when going to longer wavelengths. This is illus- 
trated by the inset in Fig.^where we show a Kurucz model 
atmosphere with Teff=7250 K before and after reddening 
with standard interstellar extinction. The extinction in the 
UV and visible is substantial but the flux level at wave- 
lengths beyond 2 /im is hardly affected. Grey extinction 
is due to dust grains much larger than the wavelength at 
which they absorb. Typically these grains should have a 
diameter of at least several /im in order to cause wave- 
length independent extinction up to 4 /im. It is unlikely 
that a large amount of such big grains grow in stellar out- 
flows. Moreover, the presence of these large dust grains 
would have strong implications for the understanding of 
the complete nebula, most notably the nebular mass bud- 
get. These implications are discussed in Sect. 18.31 

The second condition, concerning the extinction along 
the line of sight compared to the average extinction, 
is more debatable. The mid-IR images as p ublished by 
iMeixne r et all l)l997|) : iDaval et alJ l| l998^ and ' jura et"al] 
l,2000ii) show clear deviations from spherical symmetry. 
These asymmetries are due to a non-spherical distribu- 
tion of the dust around the star. We can write a simple 
expression for the flux observed from the system. 

•^obs = (1 — Vlos) -^atm + ?/avcr ^ -?'"atm (1) 
UV+visiblo IR 



I I 1 I I I I I I I I 1 I I I I 1 I 1 I I I I I 

6000 6500 7000 7500 
Te,, [K] 

Fig. 4. The ratio of the observed integrated flux over the 
integrated flux from the Kurucz model atmosphere scaled 
to the near-IR observations versus the effective temper- 
ature of the model atmosphere. The hatched area indi- 
cates the range in which this ratio equals unity within the 
uncertainties. We indicate th e position of the correspond- 
ing spectral types taken from lde Jaeer fc NieuwenhuiizerJ 
(1987) for supergiants along the top abscissa. 

•^obs _ _ 
^ — -L ^los r ^avcr : 

atm 

where rjios is the fraction of the stellar flux absorbed in 
the line of sight towards the earth and jyavcr is the frac- 
tion of the stellar flux absorbed, averaged over all lines 
of sight from the star. In this we assume that the energy 
absorbed in the CSE is re-emitted isotropically in the IR. 
The energy budget equation is further constrained by the 
observational fact that ?7avcr — 1-25 (1-ryios). This yields 

with Eq. © J^obs/^atm = 2.25-2.25?7los. 

In case ?7ios < ijavei (i-e. when we are viewing the 
star relatively unobscured) J^ohsl ^atm is larger than unity. 
This is due to light that the star emits in directions away 
from the direction of the earth. When this light is ab- 
sorbed in the CSE and partly re-emitted towards earth 
we measure extra flux. This would allow for a lower Tcff 
than the 7000 K we derive from Fig.0] However the obser- 
vations put strong constraints on the importance of this 
effect. Mathematically ryios and ?7avcr can of course be in 
the range 0—1, but in practise the UV and visible photo- 
metric data reveal a reddened star and ?7ios must be larger 
than TV (^^0.33). This is further corroborated by the fact 
that 56 per cent of the flux is emitted in the IR part of 
the spectrum. In Eq. Q a value of yyios lower than 0.33 
requires ??avor to be over 0.85. This corresponds to the un- 
likely situation in which the star is completely surrounded 
by an opaque CSE with a small hole through which we 
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happen to view the star unobscured. Reasonable values of 
7710S are within 0.44—0.6 which corresponds to J^obs/-^atm 
= 1.26-0.9 and Teff=6800-7500 K. 

It is i mportant to note that both iMeixner et al.l 
l|l997t) and IPaval et alJ ((l998l) have constructed models 
for HD 56126 in which 7710s > 7?aver in order to reproduce 
the observed mid-lR morphology. This yields a Toff in the 
high end tail of the range mentioned above. Likewise, in- 
terstellar extinction can be considered as an added contri- 
bution to ryios also yi elding a higher effective temperature 
for the central star. iBogaertI l|l994L chapter 3) has ob- 
served stars in the vicinity of HD 56126 and estimates 
the interstellar contribution to the reddening to be E(B- 
V)is==0.04. We conclude that a high effective temperature 
in the range 7000—7500 K fits best with all observed prop- 
erties of HD 56126. In the following we use T eff=7250 K 
as derived bv lVan Winckel fc Revnierl l)2000|) . 

4.2. Feature strength 

What is also directly apparent from Fig.^is the large frac- 
tion of the light emitted in the IR features. We have indi- 
cated the HAC/PAH features and the "21" and "30" /zm 
features. Together these features carry 40 per cent of the 
total IR flux (see also Table |3I). Any self-consistent radia- 
tive transfer modelling of the nebula must take absorption 
and emission by the feature carriers into account. The rea- 
son for this is twofold. First, the feature carriers cause an 
important part of the reddening of the star. Thus, when 
comparing model results to the observed photometry these 
are essential to assess the validity of the model. Second, 
in the CSE these components absorb a sizeable part of 
the stellar light and thus determine, in part, the temper- 
ature structure of the CSE as a whole. The temperature 
structure of the CSE is a very important parameter that 
translates directly into observables like the mid-IR spec- 
trum and mid-IR surface brightness. 

The strength of the IR features as seen in Fig. also 
exposes the enigma of this nebula. Several of the ob- 
served features are supposedly carried by materials com- 
posed of elements with low abundances. The key ques- 
tion is: "Is it possible that trace species can cause such 
strong features?". Note that this question applies not 
only to TiC concerning the "21" /im feature but sim- 
ilarly to the "30" /im feature which is produced by 
MgS. While Mg and S are both more abundant than 
Ti they are relatively unimpor tant in comparison to C 
l|Van Winckel fc Revnierl l2000(l that causes the contin- 
uum. 

5. Model 

The purpose of this study is to construct a detailed model 
of the CSE of HD 56126. In our modelling we focus on 
two key questions: 

— The composition of the CSE, i.e. the amount of the 
various dust components present. As we have already 



outlined in Sect, ^we observe emission from several 
components for which various carriers have been pro- 
posed. However, which of these components are present 
has not been firmly determined or the amounts of these 
carriers in the nebula have not been quantified. 

— The location of the dust. The latter is important 
since there is considerable disagreement in the liter- 
ature regarding the evidence for a superwind phase 
(with mass los s rates p ossibl y exceeding 1 0" ""^ Mf7> /vr). 
iMeixner et ail l|l997f) and iDaval et all l(l998l) de- 
scribe the mass-loss history with a low mass-loss 
phase followed b y a superwind phase. Recently, 
iKwok et alJ l)2002|) have obtained high angular res- 
olution mid-IR images and find no evidence for a 
large- scale sudden eject i on of material. On the other 
hand. fvon Helden et alJ ((200nt) have proposed that the 
"21" /im feature is carried by TiC clusters and that 
the formation of these clusters requires an short phase 
of extreme mass loss. 

Note that the observed mid-IR emission is not spher- 
ically symmetric but shows three intensity maxima 
along the limb -bright ened inner edge of the dust shell 
llMeixner et alj.1997,: .Daval et al.iri998t lJura et al.ll2000l: 
iKwok et alJl2002t) (see also Fig. EJ. We ignore this non- 
sphericity in our model because doing so greatly reduces 
the number of free parameters without compromising our 
analysis of the dust composition and density as a function 
of distance from the star. Because the nebula is largely 
optically thin, the deviations from spherical symmetry do 
not influence our conclusions. To confine the radial extent 
of the dust we compare the model results with the az- 
imuthal averaged intensity profiles derived from the mid- 
IR images. 

To model the absorption and emission of the circum- 
stellar dust we use the dust radiative transfer program 
MODUST. This code solves the radiative transfer equa- 
tion in spherical geometry subject to the constraint of ra- 
diative equilibrium. This yields the temperature of the 
dust. The code allows for the presence of several different 
dust components of various grain sizes and shapes. The 
density of the grains as a function of distance can be pre- 
scribed in various ways. MODUST yield s both the model 
SED and int ensity ru a ps. W e refer to iBouwman et alJ 
(F2OOO) and .BouwmanI <)200l[l for a description of tech- 
niques used in MODUST. The following parameters for 
modelling the HD 56126 system are used: 

— The input spectrum of the central star is g iven by 
a lin e blanketed Kurucz model atmosphere ijKurucd 
Il99l[l . We use Toff=7250 K and log(5)=1.0. With these 
parameters and the observed total integrated light of 
33.6 10~^^ W/m^ the radius and luminosity of the star 
are related to the distance by: 

i?^ = 20.5d and L* = 1051d^ (2) 

where i?* is the stellar radius in units of solar radii, L^, 
is the stellar luminosity in units of solar luminosities 
and d is the distance to the earth in kiloparsecs. 
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The location of the dust is constrained by mid-IR 
imaging. The intensity maximum which is observed at 
~1.1" from the star corresponds to the inner edge of 
the dust shell in the case of optically thin emission, 
which is certainly satisfied at IR wavelengths. Taking 
into account the point-spread-function (PSF) of the 
telescope we find that the inner edge is located at 
~1.2". Because dust emission is clearly observed in all 
mid-IR images out to a ~2.6" radius this constitutes 
a lower limit for the extent of the dust shell. However, 
the upper limit of the extent is not well constrained by 
current imaging and in the following we will treat the 
outer radius as a free parameter. 

The dust density [p) as a function of distance from the 
star (r) is specified in the form: 



(3) 



P = Po — 



where po is the dust density at the location of the inner 
edge of the shell (rp) and the power p determines the 
steepness of the density profile. For a time- independent 
mass loss p=2. If the mass-loss rate has increased over 
time p >2. 

The grain size distribution for each separate dust com- 
ponent can be specified with a power-law in the form: 



n{a) oc a (amin < a 



< 



(4) 



where a is the grain radius, n{a) it the number density 
of grains with radius a, q is the power and and 
Omax are the minimum and maximum grain radius, 
respectively. 

— We choose to keep the dust composition constant as 
a function of distance from the star. This means that 
the relative fraction of each grain material component, 
specified by chemical composition and size is fixed 
throughout the nebula. The dust shell consists of sep- 
arate grain populations, i.e. each dust grain contains 
one material, implying that separate species can have 
different temperatures at the same location. 

6. Optical properties 

In this section we discuss the optical properties of the 
various solid-state components that we use to model the 
HD 56126 system. Some of these materials have been well 
characterised in the laboratory over the broad wavelength 
range required for radiative transfer modelling. For other 
materials these data are not available and we describe how 
we treated the optical properties of these materials. 

6.1. Carbonaceous compounds 

There is general consensus that the mid-IR continuum in 
these sources is due to carbon-based dust grains, possibly 
with inclusion of hydrogen. The UIR bands at ^3.3, 6—9 
and 10—17 ^m, indicated in Fig. and |21 correspond to 
the vibrational resonances of aromatic and aliphatic C-H 



and C-C bonds. Despite this general consensus, the ex- 
act nature of the carrier is under much debate. This is 
partly due to the many different structures in which car- 
bon can appear. There are three possible bonding types 
between the carbon-atoms: sp^, the aromatic sp^ bond, as 
in graphite, and the aliphatic sp"^ bond, as in diamond. 
It is due to this diversity of bonding types of carbon and 
the resulting diversity in material properties of the vari- 
ous carbonaceous compounds that carbonaceous dust in 
space has evaded firm identification. In the following we 
consider amorphous carbon and hydrogenated amorphous 
carbon. 



6.1.1. Amorphous carbon 

The most prominent dust component that we observe is 
the mid-IR continuum. It carries ^--^60 per cent of the 
IR flux. This emission is generally attributed to amor- 
phous carbon (a-C). We use the values of the complex 
refractive index of am orphous carbon as published by 
IPreibisch et alJ l(l993l) . These data range in wavelength 
from 0.1 to 800 /im, which is sufficient to properly model 
the absorption and emission properties. 



6.1.2. Hydrogenated annorphous carbon 

The emission between 3—4 and 11 — 17 /im is indicative of 
vibrationally excited C-H bonds. This indicates that the 
CSE contains a compound of carbon and hydrogen. When 
looked at in detail both emi ssion from H-atoms bound to 
sp^ a nd sp^ sites is observed ijBuss et al.ll990tlKwok et all 
I2OOII) . Such a mixture of hydrogen bonds is generally 
found in hydrogenated amorphous carbon (HAC). HAC 
is a likely condensate in H-rich and C-rich circumstellar 
ejecta. Therefore we include HAC in our model. 

There is extensive literature on the properties of HAC 
in the material and surface science literature, where it 
is commonly referred to as a-C:H or diamond-like amor- 
phous carbon. HAC films have received much attention as 
optical coatings because they can be manufac tured to be 
relatively transparent in the 2—6 /im range; see lRobertsonI 
l)2002h . As a result of this interest, such films have b een 
studied in great detail: we refer to iRobertsonI l|2002l) for 
an excellent review of HAC properties. In an astrophys- 
ical context, the properties of HAC have been exten- 
sively discussed by Duley and co-work ers fe.g. lDulevlll984 
IScott et a,lJIT997Hcrishko etTilEool. Here we only dis- 
cuss the optical properties relevant for the current study. 
Note that because of its application as optical coating the 
material research has been focused on those materials with 
low IR absorption. The more graphite-like materials, with 
high IR absorption levels (see below) , have received much 
less attention. 

HAC is a semiconductor. Its structure can be repre- 
sented as sp^ bonded clusters surrounded by sp"^ regions. 
The optical properties from the near-UV to the near-IR 
are determined by the 7r-electrons of the graphitic sp^ 
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bonds. This is due to the fact that the 7r-electrons are 
weaker bound than the a-electrons of the sp^ hybridis- 
ation. By varying the formation conditions, the sp^/sp^ 
ratio can be varied , which in turn allow s to tune the ab- 
sorptivity of HAC ijPischler et alJll98^ . It is shown that 
the near-UV to visual absorption is determin ed by the size 
of the aromatic domains in th e material (e.g. lAngus et al.l 
Il986t lRobertsonlfl986l l2002() . HAC contains a range of 
sizes of the aromatic domains and the level of absorption 
depends on the size of the largest aromatic clusters present 
in the structure at the concentration of a few per cent. 

In the visual range, HAC absorptivity has an empiri- 
cally determined dependence on photon energy and i s usu- 
ally e xpressed in the form of the so-called Tauc law 



a{E)E = AT (E-Et) 



(5) 



where E is the photon energy of the incoming light, a{E) 
is the absorptivity as a function of energy. At is a con- 
stant that determines the slope and Et is the Tauc gap. 
Et and the optical band-gap (i.e. the energy level above 
which photons couple efficiently to the valence electrons) 
are directly linked to the average size of the sp^ domains. 

At sub-band-gap energies the electronic absorptivity 
declines exponentially and is commonly expressed in an 
Urbach law: 



a = au exp 



E 

E\j ' 



(6) 



where Eu is the Urbach energy that determines the width 
of the tail, while a\j is a constant that determines the 
level. The width of the tail, i.e. Eu, is determined by the 
disorder in the material. Hydrogenated amorphous carbon 
is a heterogeneous material with a range of bond types and 
sizes of sp^ bonded domains. As a consequence, the tail is 
broad. Typical values of the important parameters in the 
above equations are: — 300 cm~^/eV, Et — 0—2 eV 
and Eu ~ 200 meV. 

In the IR range, the contribution of the electrons to 
the absorptivity is in general small, unless one goes to 
materials with very large aromatic domains, the limiting 
case of which is the conductor graphite that exhibits a very 
strong electronic continuum in the IR. The IR absorption 
of HAC is dominated by vibrational resonances. 

We combine the observed optical properties of different 
HAC materials from various authors in order to determine 
the dependence on structural parameters from the UV to 
the IR. We use the H/C ratio as the parameter that de- 
termines the material properties. Note, that similar varia- 
tions in the sp^/sp'^ ratio as described below are observe d 
in pure amorphous carbon films (e.g. iRobertsonI |2002|) . 
However, here we study the effect of H/C because the 
astronomical spectra clearly show (through the features 
between 10 and 17 /im) that there are H-atoms bound to 
th e C-atoms. 

ICompagnini et all l)l995|) have studied the optical and 
IR absorptivity of HAC films as a function of hydrogen 
content. These authors list the parameters for both the 
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Fig. 5. The optical parameters of HAC as a function of 
H content. The panels show, from top to bottom: 1) the 
fraction of carbon bonded by sp^ bonds. 2,3,4) The inte- 
grated strength of the 11 — 17, 6—9 and the 3.3-1-3.4 /im 
bands, respectively. 5) The Urbach energy (see Eq. ©). 
6) The Tauc energy (see Eq. Q). 

Tauc law in the visible and the Urbach tail dominating 
the electronic continuum in the IR. For the HAC materi- 
als concerned here, hydrogen incorporation determines the 
structural parameters because hydrogen inclusion favours 
sp^ bonding. Hydrogen incorporation therefore limits the 
size of the graphitic domains. Thus, with increasing hydro- 
gen content the electronic absorption decreases and the 
optical gap increases. This can be seen in Fig. [SJ where 
we give the optical parameters of HAC as a function of 
hydrogen content. 

A clear correlation is found between the optical gap 
and the strength of the integrated aro matic p lus aliphatic 
CH stretching feature around 3.3 /im l|Dischler 1987) . We 
use this relation to link the strength of the CH stretching 
mo de to the electronic c ontinuum level (cf.. Fig. |SJ). 

iBounouh et alJ l)l995il discuss the influence of anneal- 
ing on the IR vibrational spectrum of HAC. Annealing re- 
moves H and as a result produces more and more graphite- 
like HAC. We use the strengths of the vibrational bands 
at 3.3, 3.4, 6—9 and 10—17 fim for their series I material 
to empirically determine the strength of the bands as a 
function of H/C ratio. 
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Fig. 6. Absorption coefficient {a) of HAC as a function of 
wavelength for changing values of hydrogen content. For 
details concerning the method used to derive the absorp- 
tivity, see text. 



HAC material also exhibits a far-IR continuum due 
to vibrations. Unfortunately, these contributions have not 
been studied well because they contain lit tle diagnostic 
value for the material properties. Following ISchutte et alJ 
lllflfl."^ . we estimate the level of the far-IR vibrational con- 
tinuum from the vibrational resonanc es observed in aro- 
matic molecules. iMoutou et al.l l)l996|) have measured the 
far-IR spectra of a variety of different PAHs. We use the 
resonances these authors find to define an average far- 
IR continuum. We apply the following procedure to con- 
struct this average. We take the wavelength and absorp- 
tion cross-section per C- atom for each of the 125 reso- 
nances listed in Tab. 3 of iMoutou et al and con- 
volve these with a Gauss function with a fuU-width-at- 
half-maximum of 80 cm"-'^. After adding the contribu- 
tions we renormalise by dividing by the total number of 
molecules (40) in the database and converting to absorp- 
tivity values in cm^^ using a typical value for the density 
of 2 g/cm^^. The far-IR contribution that we obtain in 
this way is not completely smooth but shows some struc- 
ture (see the 20—50 /im range in Fig. EJ. This spectral 
structure is most likely due to the fact that we included 
a biased sample of 'only' 40 molecules in a limited size 
range to construct the average absorption. A larger num- 
ber of molecules with a wider range of molecular struc- 
tures will probably fill in the gaps and yield a smoother 
far-IR continuum. Since this spectral structure is weak 
compared to the near- and mid-IR features and is lost 
in the emission from other dust components in the cir- 
cumstellar spectrum, we have refrained from (arbitrarily) 
smoothing it out. The optical gap and the sp^ fraction 
are correlated as explained above. The measured relation 



for a large number of HAC and a-C materials is given in 
iRobertsonI (l200l Fig. 58) fro m which we derive the sp 
fraction (cf., Fig. EJ. The final estimate for the far-IR vi- 
brational continuum is the product of the sp^ fraction and 
the average far-IR "continuum" as measured from PAHs. 
The estimate that we obtain is a lower limit to the real vi- 
brational contribution to the far-IR absorptivity because 
we only consider the contributions due to the graphitic 
(sp^) fraction in HAC, while in reality the sp'^ skeleton 
will also contribute. 

We show the dependence of the various parameters 
on the hydrogen incorporation in Fig. [3 The strong in- 
crease in the Tauc and Urbach energy with H/(H-|-C) is 
apparent. Likewise, the 3.3 and 3.4 /ini CH stretch modes 
increase with H/(H-t-C) fraction. Notice how the strength 
of the C-H out-of-plane bending modes between 11 and 
17 /im (in Fig.^lrepresented by I12) increases with decreas- 
ing H content. This seemingly contradictory behaviour is 
caused by the fact that in absolute terms the H bonded to 
sp^ sites increases during annealing. In this process some 
of the H bonded to the sp^ sites is converted into H bonded 
to sp^ sites. The absorptivity per sp^ C-H bond is much 
larger than per sp'^ C-H and the result is an increase of 
I12. In contrast the aliphatic CH stretching modes around 
3.4 /im are stronger than the aromatic modes (3.3 /im) 

The wavelength dependent absorption cross-sections 
are shown in Fig. El The main effect is that the com- 
plete continuum level from the UV to the IR decreases 
with increasing hydrogen content. Because the UV/ visible 
part of the spectrum decreases less than the IR, HAC 
with a higher H concentration will be warmer. While the 
continuum increases and the CH stretching modes at 3.3 
strengthen with increasing H content, the mid-IR reso- 
nances between 6—9 and 10—17 /im diminish in relative 
strength. In the following we refer to HAC with H/ (H-l-C) 
= X as HAC^:. 

6.2. Coal 

ICuillois et al.l l)l996(l have carried out a spectroscopic com- 
parison between the emission features of carbon-rich post- 
AGBs and terrestrial coals . This comparison looks very 
promising. Of course, the CSE of post-AGB objects do 
not contain coals as such since coals are the product of 
the biological cycle in the earth-ecosystem. However, it 
may be worthwhile to consider materials with a very dif- 
ferent formation mechanism as analogues to circumstellar 
materials. 

In stark contrast to HAC, the optical properties of coal 
have been studied much less and no spectra that cover the 
full wavelength range from the UV to the IR for one single 
sample are available. Moreover there is a large spread in 
the derived complex refractive i ndex of coals; see for ex- 
am ple the comparisons made in iManickavasagarnl l)l993f) 
and lBhattacharva et alJ l)l996^ . This large spread is prob- 
ably due to the wide variety of coal structural proper- 
ties, even for materials that are labelled the same. When 
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Fig. 7. A comparison between the optical properties of 
HAC as de rived from our mode l and "standard coal" as 
taken from iGuillois et all l|l994l) . The continuum of coal 
lies between the H/(H+C)=0.10 and the H/(H+C)=0.25 
curves. 

studying more controlled and limited samples of coal sim- 
ilar trend s as described in S ect. 16.1.21 are recovered. For 
example, llbarra et al. 1 lll99fil) find a very similar depen- 
dence of the aromatic fraction on the hydrogen content 
as found in HAC. Likewise, similar behaviour for the 
strength of the 3.3 versus the 3.4 /im band as a functio n 
of aromaticity (towards anthr acite) is found l)Cernvll996|) . 
iBrewster fc Kunitomd l)l984(l find a significant increase in 
the free electron absorptivity upon graphitisation, again 
the same as is found in HAC. Unfortunately, there are no 
coal data available that cover the UV to the IR obtained 
using one single sample. However, from the similar be- 
haviour found in coals and in HAC wc conclude that our 
HAC paramctrisation applies to coals as well. 

One important difference between coals and HAC is 
the inclusion of atom impurities, primarily O and mineral 
phases in coals. Of course, our HAC model will not be able 
to reproduce the effects of these impurities. The effects of 
impurities in the CSE of carbon stars must however be 
limited because impurities have clear spectroscopic signa- 
tures like strong bands due to 0-H (~3 /im) and C=0 
groups (5.8 ^im) l)Puente et alJll99 8V Such signatures are 
not detected in the astronomical spectrum. We conclude 
that in so far coal can be used as an analogue for astro- 
nomical carbonaceous grains, its spectral properties can 
be well represented by those of HAC. 

In Fig. 13 we show a comparison between the optical 
properties as derived from our model and t he optica l prop- 
erties of "standard coal" taken from . Guillois et al ]ri994). 
Their coal optical properties are not derived from one sin- 
gle sample and the data from 0.7 to 3 /im are not mea- 
sured but interpolated. The interpolated optical proper- 



ties in coal are higher than the Urbach tail in HACo.25- At 
the same time , the IR absorption properties adopted by 
IC^iillois et alJ l|l994l) are much weaker than the HACo.25 
features. This difference will give rise to a higher temper- 
ature for the "coal" than for HAC. 

6.3. Magnesium sulfide 

iBeeemann et al.l l|l994|) have published the measured val- 
ues of the complex refractive index [n and k values) in the 
10-500 /im range of Mg^Fe(i_:r)S with a;=0.0, 0.1, 0.5, 
0.75, 0.9. They find little difference between the samples 
with a;=0.75 and a;=0.9 and these authors suggest that 
the a;=0.9 sample is a good substitute for MgS. In order 
to reliably model the MgS contribution we need the optical 
properties in the UV to near-IR range as well. These data 
are unfortunately not available. We have opted to assume 
a grain size (amodci) and calculate the mid-IR absorption 
cross section from the n and k values in a continuous dis- 
tribution of ellipsoids (CP E) shape distribution following 
iBohren fc HuffmanI l|l983[ Chapter 12). We use the CDE 
distribution because this yields an emission profile clos- 
est to th e observed profile in many carbon-rich evo lved 
stars (e.g.lBeeemann et al.ll994HSzczerba et al.ll999| ). see 
iHonv et alJ ll2002^ for an elaborate discussion. By dividing 
the obtained cross section by the geometrical cross section 
we obtain the dimcnsionless absorption efficiency Q(A). In 
the wavelength range to 1 /im we assume Q=l, from 1 
to 2 /im we decrease Q linearly to and from 2 to 10 /tm 
we take Q=0. 

The assumptions above are of course unrealistic and 
will (most likely) overestimate the visual absorption for 
a given grain size. However, this docs not influence the 
derived mass as we can easily derive t he MgS grain tem - 
perature from the emission profile (Seelionieriniooi). 
In thermal balance the MgS temperature is determined 
by the ratio of visual over the IR absorption efficiency. 
In general the IR absorption efficiency decreases more 
rapidly with decreasing grain size than the UV and visible 
absorption efficiency. Therefore, using the MgS tempera- 
ture derived from the peak of the feature, there exists a 
grain radius (acSE) for which the absorption equals the 
observed emission in the MgS band. This value of acsE 
will be smaller or equal to the initial assumed grain size 
Omodoi- The mass that we derive does not depend on the 
grain size since the individual grains are optically thin at 
30 /im and therefore the emission scales with the mass 
(see also Sect. I7.4|l . 

6.4. Titanium carbide nano-crystals 

The identification of TiC as the carrier of the "21" /tm 
feature is based on the wavelength coincidence of a strong 
resonan ce found in TiC n ano-crystals with the "21" /im 
feature llvon Helden et al ](2000). The optical properties 
of such nano-crystals have so far not been determined in 
the laboratory. The optical properties of TiC in the bulk 
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have been studied extensively jLve & Loeothetis 


196(1 ; 


Alward et al.)|l975l Ihara et al.l|l975l ll 


jvnch et al. 


198C; 


Pfliiecr et al.' 'l984- 'Koide et al.' 'l990l 


19931: iDelin et all 


1996i:.EsDinosa-Maeaha et al..2002t). Althoueh these stud- 



ies are focused on the high energy range (the longest 
wavelength presented is 12 /im), they consistently find a 
large free ele ctron contribution to the opt ical properties in 
the infrared. iHenning &: Mutschkd l|200 present the first 
mid-IR spectroscopy of bulk TiC. Their spectrum indeed 
shows a strong mid-IR continuum due to free electrons. 
They find no resonance at 20.1 /xm and only a very weak 
resonance near 19 /im, much weaker than the resonance 
found in nano-TiC. 

In principle one can predict how some changes in 
spectral behaviour occur when comparing bulk optical 

J iroperties with those o f small (<200 atoms) clusters 
Kreibig fc VollmetjliggRl) . One of the effects is the limita- 
tion on the electron mean free-path l|Kreibig Sz FragsteirJ 
This effect is caused by "free" electrons that are 
scattered off the surface of the nano-crystal with a fre- 
quency that exceeds the frequency of collisions in the 
bulk material. In order to quantitatively calculate changes 
in the optical constants with respect to the bulk opti- 
cal properties requires separating the contribution of non- 
conducting electrons to the optical absorption from the 
free-electron contribution. For example, a detailed descrip- 
tion of how to apply this correction for silver is given in 
I Alvarez et al.l |l997.1 . This separation is based on fitting 
the optical properties over a wavelength range where the 
free-electron contribution to the absorption provides the 
dominant source of opacity. However, in TiC there is no 
wa velength range wer e this is clearly the case. For exam- 
ple |Koid^t^y(f^93^ find (besides free-electron opacity) 
inter-band transitions, due to non-conducting electrons, 
up to 12 /im; the longest wavelength they measured. As a 
consequence, the estimated contribution due to free elec- 
trons derived in one wavelength range yields results that 
clearly over- or underestimate the free-electron contribu- 
tion in other wavelength ranges. 

Another observed effect is the strengthening of the 
20 /im resonance. If we translate the strength of the 20 /im 
feature in nano-crystalline TiC to bulk samples this res- 
onance would stand out clea rly in the spectrum. Instead , 
a feature is barely detected ijHenning fc MutschkjIioOll) . 
This immediately shows that bulk TiC cannot explain the 
observed "21" /im feature. Therefore, if TiC is the car- 
rier of the "21 " /im feature the TiC must be in grains 
or in poly crystalline domains smaller than a few hundred 
atoms. 

Other quantum effects may also be important in nano- 
crystalline materials. The localisation of the electrons may 
result in a material that behaves like a semi-conductor or 
insulator in the nano-crystalline phase even if the bulk 
material is a good conductor. S uch behaviour has been 
observed for example in mercury i Busani et al.lll99^ . 

The issues outlined above leave it impossible to derive 
reliable optical constants for TiC nano-crystals in a simple 
semi-empirical way. It is clear that these issues can only 



be satisfactorily resolved by laboratory measurements of 
the optical properties of nano-TiC. Still, for our radiative 
transfer calculations some assumptions about the optical 
properties of TiC must be made. For the emission of the 
TiC grains we have adopted an approach similar to the 
one outlined above for MgS. We use the absorption pro- 
file and strength for Ti C nano-crystals as presented by 
Ivon Helden et al.l l)2000(l and assume a constant level for 
the UV/ visible absorption efficiency. This yields a TiC 
mass as a function of the UV/ visible absorption level. 



7. Radiative transfer modelling 

7.1. Initial model 

We first model only the dust continuum due to amorphous 
carbon. We use the parameters as described above. As a 
starting point for the fitting procedure we use a constant 
mass loss and outflow velocity profile (p cx r"^) during a 
period of 10000 years. This corresponds to an envelope 
with radius of ^^10". We use a single grain size of 0.1 /xm. 
Amorphous carbon grains of this size reach the highest 
temperature in radiative equilibrium, given the radiation 
field of the central star. We vary the density at the inner 
edge of the envelope (po) to fit the emission at the short- 
est wavelengths of the SWS spectrum. The emission at 
the shortest wavelengths is sensitive only to the warmest 
dust which allows us to derive the mass of the warmest 
amorphous carbon grains. We choose 9—10 and 18 /im as 
reference points. Below 9 /im the data are noisy and/or 
the contributions from the star and the features are more 
difhcult to estimate while at longer wavelengths the emis- 
sion is partly due to cooler grains. 

7.2. Cold component /extent of the nebula 

The model that we construct in this way over predicts 
the flux levels beyond 40 /xm by ^^10— 20 per cent, i.e. 
the model contains too much cold material. At the same 
time the modelled 11.9 /im and 18.2 /im spatial inten- 
sity profiles drop more steeply than observed within the 
closest 2.6 ". The same has been note d at shorte r wave - 
lengths bv iMeixner et all l)l997|) and lJura et all l)2000|) . 
iMeixner et alJ l|l997|) use p oc r~^ (p=l) in order to fit 
their mid-IR images. However, this introduces even more 
cold dust located far from the star. We compare the ob- 
served 18.2 /im intensity profile and the model profiles 
of the p=l,2 models in Fig. ^1 below. The p—2 model 
peaks too much to the centre of the nebula and produces 
too little flux in the outer nebula, while the p=l model 
is much closer to the observed intensity profile. Because 
the 'p=2'-model over predicts the flux levels at the longest 
wavelengths together with the fact that the observed in- 
tensity profllc from the inner envelope (1—2.6") requires 
p—1^ prompts us to limit the extent of the nebula. 

From an observational point of view there is no evi- 
dence that the envelope extends beyond ^^2.6". All avail- 
able images show only emission from this inner region. 
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Note, that this is not conclusive evidence because the 
available mid-IR images are all at wavelength shorter than 
25 /xm, i.e. at the Wien side of the dust emission (see 
Fig-D- This emission is extremely sensitive to the dust 
temperature and the cool dust far away from the star does 
not contribute much. Deep op tical images that show scat- 
tered light from the nebula ljUeta et alJ l2000|) show no 
evidence for emission outside the ~3" radius. Noteworthy 
in this respect is that the very high qu ality continuum irn - 
ages at 10 and 18 fiui, as presented bv iKwok et al 
are very similar in morphology and size while a dust en- 
velope that continues beyond 3" would cause the 18 fim 
images to be ~25 per cent larger than the 10 ^m image, 
where the size is measured at 10 per cent of the peak in- 
tensity. This 'growing' with increasing wavelength is not 
observed. 

We construct a model with the same parameters as 
above but with p=l. This model fits the observed 11.9 
and 18.2 /im intensity profiles much better than a law. 
From this we derive an upper limit on the outer radius of 
the envelope of 4". A larger envelope produces too much 
far-IR flux. 

Up to now we only considered 0.1 fim amorphous car- 
bon grai ns. In reality grains of different sizes are expected 
to form ("Dominik et al."l989'). We use a standard MRN 
size distribution ( Mathis et al. 197% : n(a) cx a~^-^. We 
use 0.01 /im and 1.0 fim as the smallest and largest grain 
radius, respectively. 

As described above the outer radius of the envelope 
is determined spectroscopically from the far-IR flux lev- 
els, i.e. the relative amounts of warm and cool a-C. A 
grain-size distribution affects these relative amounts be- 
cause grains of different size reach different temperatures. 
Therefore, by using a grain-size distribution the outer ra- 
dius can be further constrained. With the grain-size dis- 
tribution as described we find that with a size of 4" the 
predicted far-IR fluxes are 50 per cent too high. We flnd 
that the best fitting outer radius of the envelope is 2.6". 
We use this outer radius in the following analysis. 

7.3. HAC temperature 

We include HAC as an additional component to a-C. 
We vary H/C and the amount of HAC to best fit the 
observations. We do not aim to fit the features in de- 
tail but aim to determine if HAC reaches the temper- 
ature needed to explain the relative band strength of 
the broad emission features. With increasing H/C both 
the UV/visible and the IR absorptivity decrease. Because 
the IR absorptivity decreases more steeply, grains with 
a high H/C will be warmest. On the other hand, the 
amount of HAC required to explain the energy emit- 
ted in the features increases with H/C content because 
HAC with high H/C absorbs the stellar light less effi- 
ciently. In Table 0] we give the temperature that HAC 
grains reach at the inside of the envelope and the 
strength of the features as a function of hydrogen con- 



Table 4. HAC temperature and band strength as a func- 
tion of H/(H-|-C). H/C is the relative number of H-atoms 
to C-atoms. T^ax is the maximum temperature of the 
grains. Fg-ir^HAc/Fe-iT.obs is ratio of the integrated flux 
emitted in the plateau features in the model over the ob- 
served integrated flux in the plateau features, assuming 
equal mass of HACx and amorphous carbon. 





H/C 


Tniax 

[K] 


F6-17,HAc/F6-17,obs 


a-C 





170 




HACo.io 


0.11 


188 


1.4 


HACo.20 


0.25 


190 


1.7 


HACo.30 


0.43 


195 


1.4 


HACo.35 


0.54 


205 


0.8 


HACo.4o 


0.67 


205 


0.5 


HACo.45 


0.82 


205 


0.2 



tent. All HAC grains are warmer than a-C grains. In the 
H/(H-|-C)=0.35— 0.45 range HAC grain reach a maximum 
temperature of ^205 K. Fg-ir^HAc/Fe-iT.obs gives the in- 
tegrated strength of the plateau features relative to the 
observed strength assuming that the total HAC mass is 
equal to the a-C mass. Note, that we list the total flux 
emitted in the two plateau features. In all cases the flux 
emitted in the 10—17 fj,m band dominates and the 6—9 /im 
band is too weak compared to the observations, i.e. the 
temperature is too low. 

As can be seen from Table 01 HACo.35 is optimal in 
the sense that it reaches a high temperature while requir- 
ing 'only' equal mass as contained in a-C. To get the band 
strength of the 6—9 and 10—17 fim features right we would 
have to decrease the inner radius of the envelope. The 
band strengths are correct only when we would decrease 
the inner radius to ^0.5". At 0.5" the HAC grains reach 
a temperature of ^^290 K. However, such a small inner ra- 
dius of the dust envelope is inconsistent with the available 
mid-IR images. 

The comparison with the observations becomes worse 
in the 3 /im region. In our model HACo.35 produces no 
discernible aromatic 3.28 um or aliphatic 3.42 fim feature 
while both arc observed |Kwok ct al. 1990). We conclude 
that HAC grains in radiative equilibrium do not reach 
the temperature required to explain the observed strength 
of the 3 and 6—9 /im features relative to the 10—17 /im 
plateau features. 

7.4. Magnesium sulfide 

Before we include MgS in our model we derive the mass 
of the MgS component in a simple, relatively model inde- 
pendent, way. In the limit of grains much smaller than the 
wavelength (Rayleigh limit), the absorptivity of a grain 
scales linearly with the mass of the grain. In the optically 
thin case the observed flux can thus be expressed in terms 
of distance from the earth (d) , grain temperature (T) and 
total mass (M) alone. This can easily be derived in the fol- 
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lowing way. Consider the flux density emitted by a single 
grain (/i) of radius a and temperature T: 

fi{X,a,T{a)) = 47ra2Q(A,a)7rB(A,r(a)) (7) 

^ 3V{a)^^^7TB{X,T{a)), 
a 

where A is the wavelength, Q{X,a) is the absorption effi- 
ciency; thus 47ra^(5(A,a) is the effective radiating surface 
of the grain, ■nB{X,T{a)) is the black body flux density 
emitted per unit area and V{a) is the volume of the grain. 
In the Rayleigh limit is independent of a and there- 

fore T is independent of a and we can write: 

A(A,a,r) = 3^C(A)^B(A,r) ^ (8) 
P 

/tot(A,T) = V3^^C(A)7ri?(A,T) 
i P 

= 3^c(A)7rB(A,T), 
P 

where M{a) is the mass of the grain, p is the density of the 
material C(A) = Q^'^"-') and /tot is the total flux density 
emitted by all grains, the summation is done over all grains 
(i) and Mtot is the total mass of the grains. Integrating 
over wavelength and correcting for the distance we finally 
obtain: 

^^0 1 = / 3^CiX)TTB{X,T)dX and (9) 

Mtot = ^ Fx„_, d^p(^J^' C'(A) B(A, T) d>}j , 

where _F\o i is the observed flux integrated between A=Ao 
and A=Ai . The last term in Eq. H1U|I depends only on T for 
a given material. The peak position of the "30" u m feature 
is a good indicator of T for MgS l)Honv et al.ll2002,) and 
we derive TMgS —150 K from the observed peak position 
in the spectrum of HD 56126. Using the measured values 
p = 3g/cm^, ^^10.50 = 2.7 lO^^^ w/m^ and (i=24GG pc we 
find MMgs = 9 10-6 Mq. 

If we assume that all S is condensed in the form of MgS 
we can calculate the total corres p ondin g envelope mass 
(Monv)- IVan Winckel fc Revniersl l|2000l) have measured 
the atmospheric abundance by number of both Mg/H 
and S/H to be -410-*^, which yields a Menv of 0.06 Mq. 
The mass of the envelope is variously esti mated by fitting 
the S EP and mi d-IR imaging at 0. 54 Mq l|Meixner et al.1 
Il997ll . 0.26 Mq jPaval et al.lll998|) and 0.4 Mq that we 
derive below. We conclude that the abundances of Mg and 
S are consistent with the strength of the "30" fim feature. 

Of course, these values of AfMgS and Monv are lower 
limits. First, because we have assumed that all the MgS is 
at 150 K. In reality there will be a range of temperatures. 
The colder grains are more "hidden" and this increases 
MugS, see Eq. (|10|l . However, while such a temperature is 
too low to be compatible with the observed peak position 
of the "30" iJ,m feature, decreasing TMgS to 100 K increases 




UV/vis. -absorption coefficient [cm"'] 



Fig. 8. The mass of the TiC component as a function 
of the UV/visible absorption coefhcient (a). We show the 
TiC mass derived from the strength of the "21" /xm feature 
for different values of the UV/visible absorption coefficient 
(black line). A higher absorption coefficient requires less 
TiC. The grey diamond indicates the measured absorption 
coefficient of bulk TiC. The right abscissa shows the corre- 
sponding total envelope mass. The hatched area indicates 
reasonable values for the total envelope mass. The absorp- 
tion coefficient of nano-TiC needs to be ~20 times higher 
that bulk to be compatible with the observed strength of 
the "21" pm feature. 



MugS only by a factor 5, which still yields a reasonable 
value for M^nv Second, because we assume that all Mg- 
and S-atoms are condensed into MgS to derive Menv from 
-^MgS- We point out that there are no other Mg or S 
bearing components identified in the nebula. We conclude 
that from a simple analysis the strength of the MgS feature 
as observed in HD 56126 is compatible with the available 
Mg- and S-atoms. 

We include MgS in the radiative transfer model us- 
ing the optical properties as described in Sect. 16.31 The 
only free parameters are the MgS grain size and the 
relative mass fraction. We find that a grain size of 
augS —0.01—0.02 pm fits best with the observed profile 
of the "30" pm feature. Larger grains are too cold and 
cause the "30" pm. feature to peak at longer wavelength 
than observed. The mass of the MgS component is 4 per 
cent of the a-C, or 1.410~^Mo at d=2.4 kpc; this is in 
good agreement with the 9 10"^ M0 we derive from the 
simple considerations above, assuming TMgs=150 K. 

7.5. Titanium carbide 

As discussed above (Sect. only TiC nano-crystals are 
expected to show a strong "21" pm feature. Such nano- 
crystals will undergo a temperature spike upon absorp- 
tion of a single UV/visible photon. Coohng then occurs 
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Wavelength [|j.m] 

Fig. 9. A comparison between the observed SED (triangles) and our best fitting model (thick grey line). The model 
reproduces the dust continuum and the "30" /im feature very well. The model does not reproduce the strength of the 
6—9 /im feature and the extinction in the UV. Underneath the IR emission spectrum we give an indication of the 
contributions of the various dust components offset by -5 lO"^'^ W /ir? . 



through emission in the "21" /im feature. The mass of 
TiC is directly determined from the energy balance: the 
IR flux is equal to the absorbed UV/visible flux, which is 
given by our radiative transfer calculation. In Fig. |S1 we 
show the mass of TiC in the envelope that we derive as 
a function of the UV/visible ab sorption coefficient usin, 
the T i/H abundance of 1.3 10~* I Van Winckel fc R.evnieri 
I2OOO II. We also show the measured absorption coefficient 
of bulk TiC. It is clear from Fig.|Hlthat bulk TiC can not 
explain the strength of the feature even if all the energy 
absorbed would be emitted solely in the "21" /im reso- 
nance. In order to be compatible with the observed band 
strength nano-TiC needs to absorb 20 times more effi- 
ciently per Ti-atom in the 0.2—0.7 /im wavelength range 
than bulk TiC. This corresponds to a visual absorption 
efficiency (Qvis) of ^^8 for a lOxlOxlO-atom nano-crystal. 
Such a high value for the visual absorption coefficient is 
only possible if, due to the nanometer size scale, a strong 
electronic resonance is present in nano-TiC. This problem 
is already evident in the complete spectrum as presented 
in Fig. n] Either the carrier of the "21" /im feature is as 
abundant as MgS or it has a much higher UV/visible ab- 
sorptivity. 



8. Results and discussion 



8.1. Best fit model 



In Fig. |51 we compare the predicted and observed SED. 
Our best model fits the observations in many details. Note 
that the ordinate scale is linear while previously published 
SED fits are compared on a logarithmic basis. The peak 
of the SED, the shape and strength of the "30" fjum fea- 
ture and the shape of the "21" /im feature are very well 
explained by the model. There are two places where the 
model fails to reproduce the details of the observations: 
?) The UV flux levels are over predicted. It is clear that 
an important UV absorption component is missing in our 
model. Laboratory studies of the UV absorption proper- 
ties of MgS and TiC are therefore needed, ii) The strength 
of the 6—9 /im band is too small in the model. The dis- 
crepancy at 6—9 /xm is discussed in Sect. 18.41 

In Fig. ^] we show the comparison between the ob- 
served and predicted intensity proflles. The strength and 
position of the intensity maximum and width of the neb- 
ula are reproduced very satisfactorily by the p—\ model. 
In the outer parts of the nebula this model under predicts 
the intensities slightly. This may be due to deviations of 
the true PSF from the assumed Gaussian PSF, some weak 
background in the observation or to weak nebular emis- 
sion that we do not include in the model. The best fitting 
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Table 5. The best fitting model: parameters and derived 
quantities. 
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Fig. 10. A comparison between the observed intensity 
profile (diamonds) as a function of distance from the cen- 
tre of the nebula, the intensity profile from the best fit- 
ting radiative transfer model (line) and the profile pre- 
dicted the p=2 model (dashed line). The observed profile 
is obtained b y azimuthal averagi ng of the 18 /im image as 
published bv iKwok et alJ l|2002l) . The model profiles are 
obtained by averaging the intensity profile of the model 
shown in Fig. |^ and the p=2 model after convolution 
with a 2D-Gaussian PSF with the same fuU-width-at-half- 
maximum as the 18 /xm observation. 



p—2 model produces an intensity profile that is too much 
peaked towards the inner edge of the dust shell. 

Towards the centre of the circumstellar shell there is 
a rise in the observed intensity, which is present in all 
available mid-IR images. At 18.2 /im, this central region 
contains about 3 Jy which is ~100 times more than the 
photospheric contribution of the central star at 18.2 fim. 
Therefore, this rise must be due to the departure from 
spherical symmetry of the nebula, causing more dust to 
be present along the line of sight towards the central star 
or it is due to an additional dust component which is not 
present in our model. The nature and origin of such an 
additional dust component is at present unknown. 

We note that the emission from this region, which is 
possibly located closer to the central star than the in- 
ner edge of the dust shell in our model, cannot explain 
the plateau features and the high temperatures they re- 
quire fSect. [73)l . because the flux contributions from this 
region are much too small and the plateau features are 
emitted from a much m ore extended region, e .g. Fig. Hand 
the i mages presented in'Meixner et al.Vl997|): llDaval et alJ 
l(l998,) ; Jura et al. (2000) ; Kwok et al, (,,200^ 

We also note that our model fails to reproduce the ob- 
served intensity profiles at shorter wavelengths (10.3, 11.7 
and 12.5 /im). The observed profiles are systematically 
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broader than predicted. This is an indirect but clear indi- 
cation that the carriers of the plateau emission features at 
those wavelengths are not in radiative equilibrium, but are 
transiently heated to temperatures higher than expected 
on the basis of our equilibrium model. 

In Table El we summarise the main parameters of the 
best fit model. We find that the observations are best 
explained with a compact dust shell that extents from 
1.2—2.6". The mass of the dust contained in the shell 
is ^7 lO"* M(T, . Assuming a constant outflow velocity of 
10.7 km/s l|Knapp et alllQQSfl we calculate that the mass- 
loss commenced 2800 yr ago and halted 1300 yr ago. Using 
a value of 220, as is often used for carbon-rich circumstel- 
lar matter, for Mgas/Mjust we find that the time aver- 
aged mass-loss rate is '--^lO"'* Mq/yi. As argued below, a 
value of 600 for Mgas/Mdust is more appropriate for this 
low metallicity object, corresponding to a mass-loss rate 
of -310^'' Mq. We find that 2 per cent of the dust (by 
mass) is contained in MgS, consistent with the observed 
photospheric abundance of Mg and S. 



8.2. Previous dust models 

The circumstellar environment of HD 56126 has been 
studied before usi ng both mid-IR spectroscopic and imag- 
ing observations llMeixner et alj[l997t iDaval et alJll998t 
iHrivnak et al]l200(]Hjura et all bOOofT It is useful to com- 
pare the parameters we derive to those found before. One 
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main difference between our model and previous work is 
the fact that we treat the dust composition in consider- 
able detail. However, there are other significant differences 
as well. We find that the CSE of HD 56126 is well repre- 
sented by a compact shell of 2.6" radius. We have argued 
above that a larger outer radius will cause an increase of 
the size of the nebula with increasing wavelength, con- 
trary to the observations. Of course, a preceding phase 
with a much l ower mass-los s rate would go undetected. 
The model of ISkinner et al.l l|l995l) contains a more ex- 
tended component. The absence of this more extended 
component in our model is reflected in relatively low en- 
velope mass of 0.16 Mq that we derive. This is because 
cold material far away does not contribute to the mid-IR 
fluxes but does contain mass. This "hidden" component 
is therefore important in determining the total envelope 
mass. Imaging at sub-mm wavelengths is needed to better 
constrain the cold envelope mass. 

ISkinner et all l)l995l) find Mc„v=0.54 Mq at a distance 
of 3 kpc, with two-thirds of the mass in the extended com- 
ponent. The mass of the inner component that they derive 
corresponds well to the mass of the a-C g rains t hat wc de- 
rive after correcting for the distance. ,iDaval et a l. ( 1998) 
derive an envelope mass of 0.26 Mq at 2.7 kpc, which 
corresponds to 0.2 M© at 2.4 kpc. These authors use a 
compact model with all dust contained within 2.2" but 
with two distinct fixed temperatures. In view of the dif- 
ference in approach, the envelope mass they find and our 
mass determination are compatible. Hrivnak ct al. ( 2000) 
derive a mass- loss rate of 3.6 10~^ Mq/jt at the distance 
of 2 kpc, which corresponds to 5 10~^ M,^/vr a t 2.5 kpc, 
about half of the value that we find. lJura et~a l. (2000) 
model the density in the envelope as the result of a "dy- 
ing" wind, i.e. a decreasing mass-loss rate with time. This 
is qualitatively consistent with the po (r/Rin)~^ density 
prescription that we use. These authors find a maximum 
mass-loss rate of ~3 10~^ Mq/yi. However, they model 
the emission of the dust with a simplified amorphous car- 
bon opacity function. This yields a mass estimate com- 
patible with the warmest a-C grains (a=0.1 /zm) in our 
model. Correcting for the mass of the other grains, the 
different distance and the Mgas/Mdust ratio they use, they 
find 1.1 10""* M0/yr compatible with the mass-loss rate 
that we derived. 

8.3. Envelope mass; Mgas/^dust 

The total envelope mass that we and other authors de- 
rive scales with the gas-to-dust mass ratio. For carbon- 
rich AGB and post-AG B stars values between 200 and 
250 are often used (e.g. I.Turalll986l: iMeixner et al]ll997t) 
and in the above derivations we have used 220. However 
this value is very uncertain. Because the bulk of the dust 
is carbon-based one can expect Mgas/Mdust to strongly 
depend on the C/ ratio. The photospheri c abun dances 
determination of IVan Winckel & Revn iers^ ('2000) yield 
[C/H]=8.53-8.77 and C/O=0.7-1.4. Taking values in the 



high end of this range and assuming that all carbon is 
locked in either CO (70 per cent) or dust (30 per cent) 
we find Mgas/Mdust, c ^600, in which Mdust.c is the mass 
of the carbon-based dust component. This high value 
of Mgas/Mdust follows directly from the elemental abun- 
dance measurements. We find Mdust.c=7.4 10~^ Mq (see 
Table E)) which corresponds to Mgas=0.44 Mq. When we 
consider that the star had a ZAMS-mass of ^1.1 Mq 
and the core-mass is roughly 0.6 Mq we conclude that 
nearly the entire envelope was lost during a high mass 
loss phase of ~1500 yr. We also conclude that only a small 
amount of the dust can be "hidden" in cold grains, either 
far away from the star or in larger grains. This argues 
strongly against the importance of grey extinction due 
to large grains. Using the envelope mass of 0.44 Mq we 
find a high time-averaged mass-loss rate of a '^3 10~^ Mq. 
We conclude that HD 56126 has experienced a short pe- 
riod (1500 year) of very strong mass loss (3 10~'^ MQ/yr) 
during which the entire envelope was lost. If we adopt 
lower C/H and C/0 values this conclusion would only be 
strengthened. 

8.4. The temperature of the HAC 

One of the main discussion points concerning the relation 
between HAC (or any other solid carrier) and the 3 to 
17 ^m. features found in C-rich environments is whether a 
solid particle in radiative equilibrium is able to explain the 
observed spectra. In particular the large feature to contin- 
uum ratios found in the astronomical spectra and the rela- 
tively short wavelengths (3 /im) at which emission is found 
appear difficult to obtain in radiative equilibrium. Indeed, 
one can immediately see from Fig. Q] precisely that diffi- 
culty; the peak of the SED lies beyond 20 /im while strong 
features on a weak continuum are present in the 5—17 /xm 
range. The fundamental question is whether solid grains 
can reach the temperature required to emit efficiently at 
these wavelengths. 

We have tested this by building a detailed model of 
the optical properties of HAC. We find that in radiative 
equilibrium HAC does not reach a high enough tempera- 
ture at the location of the dust. We have tested thi s too 
for "co al" with the opacities as published bv lGuillois et alJ 
l)l994|) , even though the main difference with HAC is due 
to the interpolation these authors applied. We find that 
their "coal" reaches a maximum temperature of ~350 K at 
the inner edge of the dust shell. Which is not high enough 
to explain the relative band-strength of the plateau fea- 
tures (-500-600 K). 

This implies that the carriers that give rise to the 
plateau features are not in radiative equilibrium with 
the stellar radiation field. Most likely these features are 
due to very small dust grains or large molecules that are 
transiently heated to high temperature. From the PAH 
features that are present in the spectrum (especially at 
3.3 fj,m) it is clear that a family of aromatic molecules is 
present. The plateau features may well be carried by an 
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"extended PAH family" of aromatic molecules and clus- 
ters with attached aliphatic sides groups. 

8.5. Magnesium sulfide 

We find that the amount of MgS required to explain the 
strength of the "30" /xm feature is consistent with the 
measured abundances of Mg and S. Taking into account 
the envelope mass that we derived above (0.44 Mq for 
Mgas/Mdust=600) we conclude that only 25 per cent of 
the Mg- and S-atoms need to be condensed into MgS to 
explain the strength of the feature. Because HD 56126 
is metal-poor this implies that in other sources that are 
generally more metal-rich the strength of the "30" /im 
feature is also consistent with the MgS identification. 

8.6. Titanium carbide 

We have put TiC in our model. The optical properties 
of nano-crystalline TiC are too poorly understood to ac- 
cept or reject the TiC identification for the "21" /itm fea- 
ture. We do note clearly that the bulk optical properties of 
TiC are incompatible with the peak position, feature-over- 
continuum ratio and strength of the "21" /xm feature. The 
opacities of nano-TiC need to be ^^20 times higher in the 
visible than measured in the bulk in order to explain the 
strength of the feature. It is interesting to note that a very 
similar resonance has been found in other nano-crystalline 
metal carbides, like VC. N bC l|von Helden et alJ l200lt 
Ivan Heiinsbergen et al Of course, these species will 

not contribute much to the "21" /xm feature because of the 
much lower elemental abundance of Nb and V. Better un- 
derstanding the physics that gives rise to this resonance in 
nano-carbides may give further insight into the identifica- 
tion of the "21" ^m feature. We conclude that the question 
of the identification of the "21" /xm feature is still open 
and that nano-crystalline metal carbides remain very in- 
teresting candidates. Given the abundance of Fe and Ni, 
cubic iron-, nickel carbides (Haxonite) may be good first 
candidates for future studies. 

9. Summary and conclusions 

— We have presented a detailed study of the circumstellar 
envelope of the post-AGB star HD 56126. We present 
a large body of data from the literature, ISO spec- 
troscopy and IR imaging. From a simple energy bal- 
ance consideration we find a high effective temperature 
7250±250 K for the central star consistent with the 
temperature derived from detailed abundance analy- 
ses. 

— We have built a detailed dust radiative transfer model 
of the circumstellar envelope of HD 56126. To model 
the emission of the dust we use amorphous carbon (a- 
C), hydrogenated amorphous carbon (HAC), magne- 
sium sulfide (MgS) and titanium carbide (TiC). We 
present a detailed paramctrisation of the optical prop- 
erties of HAC as a function of H/C content. 



— The radiative transfer model explains the ISO spec- 
troscopy from 2—200 /im in great detail. The dust 
mass we derive is consistent with the dust mass as 
previously derived by other authors. We find that the 
mid-IR imaging and spectroscopy is best reproduced 
by a single dust shell from 1.2 to 2.6" radius around 
the central star. This compact shell originates from a 
short period of very high mass-loss during which the 
mass- loss rate exceeded 10~^ Mo/yr. 

— "Classical" HAC grains in radiative equilibrium with 
the local radiation field do not reach a high enough 
temperature to explain the strength of the 3.3—3.4 
and 6—9 /xm hydrocarbon features relative to the 
11 — 17 /xm hydrocarbon features. Hence we conclude 
that these features are carried by very small dust 
grains (or large molecules) which fluctuate in temper- 
ature upon absorption of a single UV/visible photon. 

— The strength of the "30" /xm feature is compatible with 
its identification with MgS. This further strengthens 
the MgS identification of the "30" /xm feature. The 
MgS temperature, derived from the observed profile 
allows an accurate determination of the MgS mass. We 
conclude that about 25 per cent of the Mg and S are 
locked up in MgS grains. 

— We find that the observed strength of the "21" /xm fea- 
ture poses a problem for the TiC identification. The 
low abundance of Ti requires very high absorption 
cross-sections in the UV and visible wavelength range 
to explain the strength of the feature. Laboratory 
studies of the UV absorption properties of TiC nano- 
crystals are therefore urgently needed. Other nano- 
crystalline metal carbides also exhibit similar reso- 
nances as found in nano-TiC and these metal carbides 
should also be investigated in the laboratory to clarify 
the nature of these resonances and the identification 
of the "21" /im feature. 
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